A singular integrable equation from short capillary-gravity waves 
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From a columnar approximation of the Euler equations of an incompressible fluid with surface 
tension, we derive in the short-wave approximation a new integrable classical 1+1 dimensional field 
theory for the motion of the surface. Together with a Lorentz invariance, this system has the novel 
feature of solutions which become multiple valued in finite time. 
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Introduction. The nonlinear and dispersive propaga- 0(1) and 
tion of surface waves in an ideal incompressible fluid 
(depth h, density a), under the action of gravity g and 
surface tension T, is a classical subject of investigation 
in mathematical physics 0, 0, Q ■ In this letter we derive 
and study a new integrable model equation from asymp- 
totic dynamics of a short capillary-gravity wave, namely 
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Here u{x, t) is the fluid velocity on the surface, x and t are 
space and time variables, subindices mean partial deriva- 
tives, 9 = (T/ah 2 g) is the dimensionless Bond number 
and v = {ST/ah) 1 / 2 . 

The dynamics of surface waves in an ideal fluid obeys 
complicated nonlinear and dispersive equations. To sim- 
plify them, multiscale asymptotic methods can be em- 
ployed. Most of the resultant asymptotic models rep- 
resent, for large t, balance between weak nonlinearities 
and linear dispersion. For instance the long-wave dy- 
namics of a low amplitude initial profile on a shallow 
dispersive fluid are well known nowadays. The models 
extend from the oldest Boussinesq systems or the ubiq- 
uitous Korteweg-deVries Q to the more recent Camassa- 
Holm equation jj] with nonlinear dispersion. In contrast, 
hardly anything is known about asymptotic models for 
nonlinear and dispersive dynamics of short-waves. For 
the most part short waves have been studied in connec- 
tion with modulation of short-wave trains |5|, Ifl I3j l3> l2| • 

In this paper we derive the model in the short-wave 
regime of surface waves, we prove that it is integrable 
and show that it leads to unusual special solutions that 
develop singular behavior in finite time. 

The short-wave limit. To define a short wave (wave 
length I, wave number k = 2tt/1) one needs to compare 
I to an underlying space scale. We use the unperturbed 
depth h as the natural reference, and thus consider h = 



where e is the parameter of the asymtotic expansions. 

Multiscale asymptotic methods are strongly based on 
the dispersion relation u)(k) and on the associated phase 
velocity v p and group velocity v g . The short-wave limit 
is meaningful if and only if those two velocities possess 
a finite limit. Then v p and v g allow to define asymptotic 
variables and to handle the nonlinear regime (lOL fll| - 

For the usual linearization of the Euler equations (with 
surface tension), the linear dispersion relation 



Tk 



w(k) = [k(g+ )tanh(fc/i)] 1/2 



(3) 



yields in the short-wave limit v p ~ (Tk/a) 1 ^ 2 — * oo. This 
not only prevents us from defining asymptotic variables 
but also infinite dispersion cannot be compensated by 
weak nonlinearities. We found that the solution to this 
problem is to employ the Green-Nagdhi conditions of lin- 
earization. 

The basic model. Green, Laws and Nagdhi [L^ . [l3L[T^| 
developed alternative reductions of the Euler equations 
leading to models having dispersion relations with good 
behavior in the short-wave limit, as demonstrated in |lfj|. 
They used three main hypothesis, namely non irrota- 
tional fluid flow, motion in vertical columns and non- 
Archimedian pressure condition. For seek of complete- 
ness, we derive hereafter the model in a simple manner 
and include surface tension. 

Let the particles of the fluid be identified in a fixed rect- 
angular Cartesian system of center O and axes (x,y,z) 
with Oz the upward vertical direction. We assume trans- 
lational symmetry in y and we will only consider a sheet 
of fluid in the xz plane. This fluid sheet is moving on 
a rigid bottom at z = and its upper free surface is 
2 = S(x,t). The continuity equation and the Newton 
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equations (in the flow domain) read 
u x + w z = 0, 

a(u t + uu x + wu z ) = -p*, 
<r(wt + uw x + ww z ) = —p* 



9°~ 



(4) 
(5) 
(6) 



where p*(x,z,t) is the pressure and (u,w) the vectorial 
velocity. 

The kinematic and dynamic boundary conditions read 



w = at z = , 
St + uS x — w — at z 

* T ' S xx 



P =Po 
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S(x,t), 
at z = S(x, t). 



(7) 
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The columnar-flow hypothesis consists in assuming that 
u does not dependent on z, hence from and Q) 

u = u{x,t) , w = —zu x . (10) 

The integration of over z from to S^x, t) then gives 

«7S(ut + uu x ) = -p x + T[(l + S£)~*]x, (11) 



p(ac,t) 



S(x,t) 



p*(x, z, t)dz — pqS{x, t). 



(12) 



Now we multiply © by z and integrate it over z to get 
S\ . , . 7.S.S,,, '/a.S- 



P" 



(13) 



Finally, elimination of p between 111 II) and (|13|) gives, to- 
gether with |JSJ and (|1(J|> . the extension of the Green- 
Nagdhi system to non-zero surface tension 



S(u t + uu x ) =- [S 3 (u xt + 



x- u 2 x )] x ~ gSS x + 
(T/a)S S xx (l + S 2 x y 3 ^ 2 ] , (14) 

J X 

S t + {uS) x = . (15) 

This constitutes our basic model. 

Asymptotic model for short capillary- gravity waves. 
In contrast with shallow water theories with dispersion 
(Boussinesq type equations), this model incorporates fi- 
nite dispersion both in the long-wave and in the short- 
wave limits. Indeed the linear dispersion relation is 

n 2 = k 2 [gh + (Th/a)k 2 ]/[l + (hk) 2 /3}. (16) 

Hence the phase velocity is bounded in the short-wave 
limit as we have 
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(17) 



This allows to define asymptotic variables [l(j 

C= (l/e)(x-vt) , T = et. (18) 



With the power series u — e 2 (uo + e 2 U2 



and S = 



h + e 2 {S a + e 2 S 2 + ■ ■ .), the basic system ([1311(113 leads to 
an equation for Uq(Cj t) which, in the laboratory variable, 
becomes our main equation ljT)l. 



Lax pair and finite-time singularities. After appropri- 
ate rescalings of the variables, one can bring equation 
into the form 



u x t = u — uu xx 



\ u l + ^u xx ul, (19) 



A being expressed in terms of the physical parameters of 
equation Q . The corresponding Lagrangian is: 



r 1 1 2 1 2 

£ = ^u x u t + -u + -uu x 
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(20) 



Equation (|19Jl is integrable with Lax pair (in usual no- 
tations): 



L = 



d_ 

dx 
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where a are the usual Pauli matrices, E the "eigenvalue" 
and 



F 2 = 1 - *2,u xx + Xu 
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xx ' 



(23) 



One of its most remarkable properties is that with F one 
builds the first non-trivial conserved quantity for all A: 



Ft 



1< ' F 



(24) 



and through the change of function from u(x, t) to 
1 



g{y,t) 



with 



VT^X 



y 



Argtanh^^-A (25) 

1 ^XX 



Fdx, 



one finds that g satisfies the sinh-Gordon equation 
1 



9yt = 



VT^X 



sinhV 1 — Xg. 



(26) 



(27) 



This is valid for A < 1 and for u xx small enough so 
that F is real. If u xx is large, a similar change leads to 
the cosh-Gordon equation, and if A > 1 one obtains the 
sine-Gordon equation. Finally, for A = 1, one obtains 
for g{y, t) — u xx / (1 — u xx ) the equation of a free field in 
light-cone coordinates. 

Whatever the value of A, it follows from the change 
of variables (x,t) ^* (y,i) that a regular g(y,t) can give 
back a singular, multivalued u{x,t) if the change from y 
back to x is not one-to-one. This happens when \g\ is 
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large enough, forcing u xx to infinity and a change of sign 
of F in the equation for y. We give an example of this in 
fig. 1, where two solutions for u are plotted from breather 
solutions of the sine-Gordon equation, one (dashed curve) 
with an amplitude just below the singularity threshold, 
so that u is still regular and single- valued, the other (solid 
curve) with an amplitude above the singularity threshold, 
which displays a swallowtail behavior. 
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FIG. 1: u(x,t — 0) corresponding to sine-Gordon breathers 
g(y, t)V\^T = -4Arctan |gf[f[fff ,c 2 + s 2 = 1, d = s/c, at 
A = 2.02, d = 0.2 (dotted curve), d = 0.35 (solid curve). 

Furthermore, this singularity and multivaluedness can 
be reached in finite time from a regular solution. In par- 
ticular, in the sine-Gordon case, one can start from a 
solution g(y, t) consisting of two breathers far enough 
from each other and each weak enough so that the change 
y —>■ x is single valued, but strong enough so that when 
they overlap, \g\ becomes large enough for the singular- 
ity to appear. In the free field case A = 1, and even in 
the sinh-Gordon case, one can replace the breathers with 
wave packets, and make them collide and give rise to this 
singularity before they disperse. In all these cases, the 
singularity in u is a swallowtail just like in figure 1. 

The (singular, multivalued) u(x,t) corresponding to 
the sine-Gordon soliton for g{y,t) is also interesting. It 
is displayed on figure 2 for the particular value A = 10/9. 
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FIG. 2: u(x,t — 0) corresponding to the sine-Gordon soliton 
g{y,t)^\ - 1 =4Arctanexp(y-M) at A = 10/9. 



For A = 0, equation ljl§|) was already discussed in 
El 13 Elll El H3, and it contains peakons. These 
peakons, which are solitons, and their scattering are qual- 
itatively very easy to study via the change of variables to 
the sinh-Gordon equation, where they correspond to sin- 
gular solutions obtained by simple analytic continuations 
of the sine-Gordon multisoliton solutions. 

For any value of A, the equation 119(1 has the symme- 
try x — > ax, t —y t/a, u — > a 2 u, for arbitrary real 
a, which is just the Lorentz group in 1+1 dimensions. 
Although the Lagrangian i|20|) does not have the appro- 
priate covariance property to give a Lorentz invariant 
action, the system, being integrable, has an infinite set 
of Lagrangians (and derived Hamiltonians and Poisson 
brackets) leading to the same equation of motion 119fl . 
and one of them, built with the invariant field u XXl leads 
to an invariant action. 

Benjamin- Feir instability. The Benjamin- Feir insta- 
bility results from resonant interaction of an initial 
monochromatic wave with side-band modes produced 
by nonlincarity. This instability, which modulates the 
carrier envelope, is well described by the nonlinear 
Schrodinger asymptotic limit. 

Following the standard approach Q we can show that 
a Stokes wave train of equation Q is unstable if 



namely any slight deformation of the plane wave expe- 
riences an exponential growth. In the case of water at 
room temperature (T = 0.074 Nm^ 1 , a = 10 3 kg m~ 3 ), 
we obtain that a short wave train is unstable for a depth 
h > 0.49 cm. 

Last but not least, the value 9 = 3/10 corresponds to 
A = 1 in ljHf|) . Precisely 

6 < 0.3 => A > 1. (29) 

Comments. The occurence of singular (even multival- 
ued) solutions in an equation derived in a hydrodynamics 
framework is interesting, especially since u is the devia- 
tion of the free surface from equilibrium. The actual ob- 
servability of the present singular behaviors would require 
a more detailed analysis of the validity of the short-wave 
approximation near the singular points which goes be- 
yond the scope of this paper. In particular the inclusion 
of viscosity, which acts strongly over small scales, will 
affect the short-waves dynamics and alter these singular 
behaviors. 

The Green-Naghdi equations can be improved system- 
atically toward higher depths |2l|. In the linear limit, 
these improvements give the higher (N, N) Pade approx- 
imants of the Euler dispersion relation @, the present 
case (11 (if corresponding to N = 1. In particular they 
always lead to a finite phase velocity in the short-wave 
limit, which grows quickly as the order of the approxi- 
mation grows to try to mimic the behavior of the exact 
formula. 
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What is remarkable is that the short-wave asymptotics 
of these improvements lead to exactly the same integrable 
equation except for different numerical coefficients in 
front of the physical quantities g, h, 9 and v. Hence we 
hope that at least some the main features of the singular 
behavior of the solutions correspond to the actual physics 
of these water waves in arbitrary depth. 

Another point is the existence of a peakon solution in 
the case A = 0. In the rescaling leading from equation 
Q to equation 119|l . this value of A is obtained only for 
9 = 1/3, where the rescaling is singular for the x, t and u 
variables themselves, so that the whole asymptotics must 
be reconsidered from the start. This value 9 = 1/3 in the 
Euler equation leads to a dispersionless system for small 
k only and it is a peculiar feature of the Green-Naghdi 
equations to be dispersionless for all k (hence for large k) 
for that value of 9. This large k feature is not inherited 
by the improvements of [2j . 

Equation Q has a Lorentz invariance, and its quantum 
version promises to exhibit new features not shared by 
the existing relativistic systems in 1+1 dimension. This 
new relativistic integrable system, with just one massive 
bosonic field and a local classical equation of motion, is 
quite intriguing. In particular, the change of variable 
which transforms it into sine- or sinh-Gordon mixes the 



space-time variable x and the field u xx and furthermore 
requires the equation of motion to be satisfied. Hence, 
it cannot be quantum mechanically equivalent to sine- or 
sinh-Gordon. For example, it is not parity invariant (par- 
ity in the laboratory frame, x + t — ► —x — t, x — t — > x — t). 
From this follows an unusual S'-matrix, in which there 
are two phase-shifts, one for the left moving particle and 
one for the right moving one. The quantum field the- 
ory and mathematical structures following from this are 
worth a detailed study for themselves in separate pub- 
lications. This could also have some analogy in general 
relativity, where space-time is, through the metric, a dy- 
namical variable and where one can go from a choice of 
space-time parametrization to another one by a change 
which can involve the metric itself. 
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